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We study the nucleation dynamics of a model solid state transformation and the criterion for 
microstructure selection using a coarse-grained molecular dynamics (MD) simulation. Our simu- 
lations show a range of microstructures depending on the depth of quench. We closely follow the 
dynamics of the solid and find that transient non-affine zones (NAZ) are created at and evolve 
with the rapidly moving transformation front. The dynamics of these plastic regions determines the 
selection of microstructure. We formulate an elastoplastic theory which couples the elastic strain to 
the non-affine deformation, and recover all the qualitative features of the MD simulation. Using this 
theory, we construct a dynamical phase diagram for microstructure selection, in addition to making 
definite testable predictions. 

PACS numbers: 



I. INTRODUCTION 

The dynamics following a quench across a solid state 
structural transition, rarely takes the solid to its equi- 
librium state Severe dynamical constraints experi- 
enced by the product inclusion within the parent crys- 
tal, determine the mode of nucleation and of subsequent 
growth. Often, solids get stuck in long-lived microstruc- 
tures, which depend on the depth of quench and cool- 
ing rateQ. For example, transformations occurring at 
high temperatures are typically accompanied by large- 
scale rearrangements of atoms; in this case the elastic- 
ity of the solidplays only a minor role in determining 
microstructure [J [3, [1]. On the other hand, at low tem- 
peratures, only local rearrangements of atoms are possi- 
ble; the resulting microstructures are largely determined 
by elasticity P, 0. These are just two of the myriad pos- 
sibilities explored by the transforming solid. Which of 
these is actually selected; in other words, can we con- 
struct a dynamical phase diagram? 

In a set of papers 0, H, H], we had explored these is- 
sues using an MD simulation Q of a model system un- 
dergoing a two dimensional square to rhombic structural 
transformation. We found that when the transforma- 
tion proceeds at a high temperature, the resulting prod- 
uct nucleus is isotropic and polycrystalline, while a low 
transformation temperature induces the formation of an 
anisotropic nucleus, roughlyelliptical, consisting of a pair 
of twin-related crystallites p. The two modes of nucle- 
ation may be denoted Ferrite and Martensite, borrowing 
terminology from the microstructure of steel m. By fol- 
lowing the nucleation dynamics in 'microscopic' detail, 
we had established that the ferrite nucleus is formed fol- 
lowing extensive rearrangements of atomic coordinates, 



large extent, preserved. This is consistent with the two 
paradigms commonly described in real materials. How- 
ever, these two limits are not mutually exclusive @; in- 
deed for intermediate temperatures, the transformation 
proceeds such that both mechanisms may operate at dif- 
ferent spatial and temporal locations [g, a feature ob- 
served in real materials[8|. Further, the different mi- 
crostructures (twinned and un-twinned) were obtained 
simply by tuning appropriate kinetic parameters @, as 
observed in the heterogeneous nucleation of colloidal 
crystals [9|]. These observations underline the need for a 
unified theory of microstructure selection describing the 
dynamics of nucleation of both the ferrite and martensite 
and the conditions in which these microstructures obtain. 

Our preliminary attempts at this unifying picture^] 
were based on the recognition (from the MD simula- 
tion) of the role played by non-elastic variables, which 
we identified with local density fluctuations. We showed 
that the coupled dynamics of density fluctuations and 
elastic strain determined the microstructure of the grow- 
ing nucleusQ. Here, we provide a more refined the- 
ory of solid-state nucleation and microstructure selec- 
tion. Based on our MD simulations, we formulate an 
elastoplastic theory which interpolates between ferrite 
and martensite as some relevant parameter is tuned. We 
discuss the dynamical origin of microstructure selection 
and exhibit a dynamical phase diagram [g^ of the final 
microstructure. 



A. Average compatibility : geometrical versus 
strain-only theories 

Geometrically, a martensite results from the require- 
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affine transformation, viz., 

R' = T(R,t)R 

where R and R' are lattice vectors in the parent and 
product structures, respectively. The constraint that the 
martensite product smoothly develops from the parent 
phase, requires that they coexist on either side of a locally 
planar interface. This translates to the requirement that 
the transformation T possesses such a plane on which 
an arbitrary vector remains untransformed[l|, [loj - the 
constraint of rank one compatibility, which is usually im- 
possible to satisfy for any single T. Nevertheless, it may 
be possible to obtain an 'approximate' planar interface 
between the parent and product structures when two or 
more degenerate variants (twins) of the transformation 
are used (Fig. 1). This interface possesses only average 
rank one compatibility [ 1 Oj so that the rotation of an in- 
terfacial vector, averaged over a coarse graining length 
along the interface vanishes. Given a T, crystallographic 
theories[10] of martensite structure list all possible inter- 
faces with average rank one compatibility (also known 
as an "invariant plane strain"). Note that at distances 
smaller than the size of the twins, the interface is not 
planar; at this scale the transformed region is in general 
not obtainable by an afhne transformation of the parent. 
Mathematically, one may define a series of interfaces with 
rank one compatibility over ever finer scales using gra- 
dient Young measures [13, [HI- The convergence of this 
series cannot be determined within geometrical theory 
alone, and needs an additional physical input. By includ- 
ing the interfacial energy cost for creating twin variants, 
the series of interfaces may be made to converge in the 
mean to a limiting interface over the Young measure. 

A parallel approach to the study of martensitic struc- 
tures was initiated by Barsch and Krumhansl[12]. In 
this framework, the martensitic structure results from 
the minimization of a non-linear, elastic free-energy func- 
tional where the components of elastic strain are used as 
order parameter (OP) distinguishing parent and product. 
Unlike geometrical theories however, this programme 
can be developed to study the dynamical evolution of 
microstructurefisl. [3, [3, The driving force for 

the nucleation dynamics of martensite from the austen- 
ite is derived from the same free-energy functional, writ- 
ten in terms of a dynamical elastic strain tensor. The 
free-energy functional supports several degenerate min- 
ima corresponding to the different variants of the prod- 
uct. Unlike the geometrical theories, these 'strain-only' 
theories impose exact elastic compatibility at all space 
and time, restricting the elastic displacements u(R, t) = 
R'(t) — R to smooth, single-valued functions. This local 
constraint automatically disallows all configurations with 
defects and regions of plasticity and assumes that the 
transformation is affine over all length and time scales. 
Integrating out the non-order parameter (NOP) strain 
components, introduces non-local interactions between 




FIG. 1: Schematic of the austenite-martensite interface : the 
(almost) planar interface (dashed line) separates the austenite 
from the martensite, which consists of two variants denoted 
by the transformations T and T'. Elastic compatibility at 
this interface is restored when the atomic coordinates implied 
by this diagram are coarse-grained over distances of order A, 
comparable to the width of the twin variants. 



The exact compatibility of strain-only theories seems 
to be at variance with the average compatibility of ge- 
ometrical theories, and indeed with experimental ev- 
idence from micrographs of the austenite-martensite 
interface fl7j. One of the consequences of our study will 
be to reconcile these two notions of compatibility. 

To recognise this as a problem, let us recall the tradi- 
tional hydrodynamic description of the long length and 
time scale dynamics of a solid. This hydrodynamic de- 
scription is written in terms of conserved variables (such 
as mass density and momentum density) and broken 
symmetry variables (such as the displacement field or 
alternately, the strain tensor) [l^, [l^. The strain-only 
theoriesflJI, go a level further, in declaring that the lo- 
cal dynamical evolution of a solid is solely via affine de- 
formations of the parent solid. The local compatibil- 
ity constraint follows naturally from this restriction on 
the space of allowed configurations explored by the solid, 
since configurations involving not only non-affine defor- 
mations, but vacancies, dislocations and other defects are 
disallowed from the dynamics. With this restriction, the 
transforming solid has no choice but to be a martensite. 
Strain-only theories are therefore incapable of describing 
the ferrite or a host of other microstructures with inter- 
mediate characteristics where non-affine transformations 
are involved Q. 

In order to provide a unified description of the dynam- 
ics of solid state transformations, we will need to break 
this impasse. One way of achieving this, is to recognise 
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structural transition. There is therefore no fundamental 
reason to restrict the dynamical variables to hydrody- 
namical variables. The coupling between the dynamics 
of these faster degrees of freedom and the conventional 
hydrodynamical variables could in principle give rise to 
a criterion for microstructure selection. The problem is 
- what constitutes the relevant faster degree of freedom. 
In the absence of an apriori fundamental principle, we 
will use our coarse-grained MD simulation to explore all 
possible deformations left out by the earlier strain-only 
theories, namely, non-affine deformations, vacancies and 
dislocations. Of these candidate fast- variables, we will 
see from our MD simulations, that the most relevant 
ones are transient and localised non-affine deformations, 
which inevitably accompany the tranfomation regions. 

An extension of strain-only theories to include the dy- 
namical effects of non-affine deformations of the parent 
lattice, as outlined in Sect. Ill, results in (i) a proper the- 
ory of microstructure selection and (ii) elastic compati- 
bility only when merged over a coarse graining distance 
of the order of the size of the twin domains. This fea- 
ture is present in our MD simulations Q . This represents 
a reconciliation between strain-only and geometric theo- 
ries. Indeed, as we show here, elastic compatibility, and 
only on the average, emerges from a dynamical theory 
unifying the nucleation of a ferrite and martensite. 



B. Elastoplastic description : microstructure 
selection and resolution of compatibility 

Non-affine or plastic deformations have been the sub- 
ject of much study in solids subject to large shear defor- 
mations. The many approaches to the study of plasticity 
include the phenomenological elastoplastic theories 21', 
[22j and the more 'fundamental' non-affine field theories 
due to Falk and Langer[23j and Lemaitre24|. Central 
to these approaches, is the decomposition of the total 
strain into elastic and plastic parts. A plastic strain de- 
velops once the local stress exceeds a yield stress. In the 
elastoplastic theories of shear deformed solids, the yield 
stress and the dynamical constitutive relations between 
stress and deformation rate (plasticity), are phenomeno- 
logically introduced or 'derived' from models of inter- 
acting dislocations [2ll. [2^ . The more recent non-afhne 
field theory approach|23l l24l| attempts at a unified de- 
scription of plasticity in crystalline and amorphous solids, 
in terms of microscopically defined shear transformation 
zones (STZ) representing local regions with high non- 
affine deformation. The coupled dynamics of STZs and 
elastic strain produces a 'first-principles' description of 
plasticity, yield and work hardening which is consistent 
with the phenomenological elastoplastic theories. 

We follow a similar program in our study of microstruc- 
ture selection in solid state transformations. Our ap- 
proach is guided by MD simulations of the nucleation 



the transformation, create local non-affine zones (NAZ) 
beyond a threshold stress. Within these zones the affine 
connection between the parent and product lattices - and 
the smoothness of the displacement field taking the par- 
ent crystal as reference - breaks down so that the product 
crystal cannot be described as a purely elastic distortion 
of the parent. One may also regard NAZs as regions of 
high dislocation density though this description is not 
particularly useful because at such high defect densities, 
the identity of individual dislocations is lost. We find that 
the dynamics of these NAZs determines the selection of 
microstructure. We then highlight four generic principles 
derived from our MD simulations, and use these to con- 
struct an elastoplastic theory for the dynamics of solid 
state transformations, in terms of a non-order parameter 
(NOP) plastic strain describing the viscoplastic nature of 
NAZ. We show that this elastoplastic theory successfully 
describes both the ferrite and martensite nucleation and 
microstructure, and display a nonequilibrium phase dia- 
gram. While describing the dynamics towards a marten- 
site microstructure, our theory reduces to a variant of 
the strain-only theory [l^, when the austenite-martensite 
interface is coarse-grained over a length scale X of the or- 
der of the typical size of the NAZs. In this way, average 
compatibility emerges from the dynamics describing the 
martensite. We follow this up with a studyf^T*], which as- 
sociates the NAZs obtained in our MD simulations, with 
distinct particle trajectories. 

The rest of the paper is organized as follows. In the 
next section, we describe our MD simulations on the 
square to rhombic transition, with special emphasis on 
the identification of non-affine zones during the solid 
state transformation. We show how local yielding is asso- 
ciated with a change in the resulting microstructure. In 
section III, we develop an elastoplastic theory of solid- 
state transformations and apply it to the specific case 
of the square to rhombic transition. We show how our 
coarse-grained theory, qualitatively reproduces the main 
features of the MD simulations. In section IV, we dis- 
cuss some implications of our study, list some unresolved 
questions and indicate directions of future work. 



II. MOLECULAR DYNAMICS SIMULATIONS 
OF A MODEL SOLID STATE 
TRANSFORMATION 

There are many real two dimensional systems which 
show structural transitions. These include confined 
molecular and colloidal sohds |2S], flux lattices ^2^, 
skyrmions in fractional quantum Hall systems [s^], mag- 
netic colloidal particles [3l| and colloids in electric 
fields (3^. Our aim is not to mimic any of these systems in 
detail, but to construct a generic model which is able to 
describe such lattice transformations in two dimensions 
using a coarse-grained potential. In addition, we would 
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the model potential. 

A simple effective model which shows transitions be- 
tween square and rhombic (a special case of the more 
general oblique) lattices, comprises of particles which in- 
teract via the potential^, 

l/2^F2(r,,) + l/6 J2 ^^3(r.,r„r,), (1) 

where is the position vector of particle i, and rtj = 
\rij\ = |rj — ri|. The anisotropic two-body potential 
is purely repulsive and short ranged, 

V2inj) = V2 i^^j {1 + acos^ 2%} (2) 

where (Tq and V2 set the units of length and energy, a is 
an 'anisotropic lock-in' parameter|^], and 9ij is the angle 
between ry and an arbitrary external axis. The short- 
ranged three-body interaction [ssj. 

V3(ri, Tj, Tfe) = V3 [fijfjk sin^ 46*^7, -f permutations] , 

(3) 

where the function = f{i~ij) = ^ '"o)^ for r^j < 
ro = 1.8(70 and otherwise and the angle 6ijk is the 
angle between the vectors and rjk- The two-body 
and three-body interactions favor rhombic and square 
ground states, respectively. Inclusion of the two-body 
anisotropic lock-in parameter a is a device to vary the 
jump in the order parameter from strongly first order 
(a = 0) to a continuous transition for W3 = 0, a ~ 1.5. 
The unit of time is (Jq\J m|v2^ where m is the particle 
mass. Using typical values, this translates to an MD 
time unit of Ips. Knowing the individual particle MD 
trajectories allows us to project time dependent atomic 
positions into time varying coarse-grained fields whose 
evolution can be monitored during the transformation. 

Both the two and three -body potentials are purely 
repulsive and therefore the system needs to be con- 
fined either in a box of fixed volume or by an external 
compression 0. In this paper, we discuss our results for 
MD simulations in the constant number, volume (and 
shape), and temperature (NVT) ensemble with periodic 
boundary conditions using a Nose-Hoover thermostat. 
We have, in addition, carried out extensive simulations 
in the constant stress (NET) ensemble with open bound- 
aries using an additional confining potential which, at 
the same time, allows for changes of overall shape of the 
crystal during the transformation. Our main results con- 
cerning the dynamics and mechanism of microstructure 
selection are the same in both the ensembles. A detailed 
comparison of MD simulation of our model system in 
various ensembles, starting from a variety of initial states 
and for the full range of the potential parameters is being 
prepared for publication elsewhere. 

The equations of motion for up to = 20000 parti- 
cles are integrated using a Verlet scheme with a time 




FIG. 2: (color-online) (a) Phase diagram in the T — V3 plane 
(with a = 0), the solid line is the phase boundary be- 
tween square and rhombic crystals. The dashed line marks 
the temperature above which anisotropic, twinned nuclei be- 
come rare and is identified as the martensite-start (or Ads) 
temperature[2| for our model. Typical product nuclei formed 
following quenches at (b) T = 0.8 («3 = 10 ^ 5.5), to obtain 
the isotropic 'ferrite' and (c) T — 0.1 {v3 = 5 —> 1.65) to 
obtain the anisotropic, twinned martensite, starting with an 
equilibrated square parent crystal composed of 12099 parti- 
cles. 

Pq — N/V — 1.1 is obtained by computing and compar- 
ing the free energies of square and rhombic lattices using 
the technique outlined in [s^ . 

We have now set the stage for a detailed study of the 
nucleation dynamics of a solid in solid, following a quench 
across the structural transition. Our effort will be to ex- 
tract general matters of principle from these simulations; 
we will highlight these as we go along. These principles 
will form the basis for the elastoplastic theory of solid 
state transformations (section III). 

A typical quench from a square to a rhombic solid 
into a region where the square lattice is metastable, initi- 
ates multiple nucleation events (at least at high temper- 
atures), making a quantitative analysis of the dynamics 
of a single critical nucleus, cumbersome. We get over this 
difficulty by introducing a nucleation seed at the center of 
the simulation box. The seeding consists of replacing the 
central particle with a particle whose size a is smaller by 
a factor S = (cto — a)/ao, we have taken 0.25 < <5 < 1, so 
as to obtain nucleation events within reasonable compu- 
tation time. Having equilibrated the seeded square crys- 
tal at large V3, we "quench" across the phase coexistence 
line by varying the coefficient of the three-body term V3 
at two different temperatures T = 0.8 and T = O.ljs^. 
While the seeding is a matter of convenience at the higher 
temperature, it is necessary at the lower temperature. 
The transformation at the lower temperature proceeds 
via heterogeneous nucleation [6|. 

In Figl2l^b),(c), we show a snapshot of the resulting 
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tively. The colors indicate the local bond-angle order 
parameter which is defined to vary from (blue) in 
the square lattice to 1 (red) in the rhombic [6|. It is 
clear from the particle position snapshots, Fig[2I^b),(c), 
that the product nucleus is isotropic for large temper- 
atures and highly anisotropic for small temperatures. 
We identify the isotropic nucleus with a ferrite and 
the anisotropic one with martensite[6\. This identifi- 
cation is reinforced by showing that the latter is twinned. 

1. Solid state transformations predominantly proceed 
via nucleation. At low temperatures, the nucleation of 
the product solid is heterogeneous and is initiated by 
'seeding' the parent. 

To follow the dynamics in quantitative detail, we com- 
pute the coarse-grainjjd local strain field using the proce- 
dure introduced in 23]. Briefiy, we compare the imme- 
diate neighborhood 17, centered around r, of any tagged 
particle (defined using a cutoff distance equal to the 
range of the potential) in the initial, reference, lattice (at 
time t — 0) with that of the same particle in the trans- 
formed lattice. We obtain the "best fit" local affine strain 
eij = Tij — Sij {dij is the Kronecker tensor) which maps 
as nearly as possible all the particles n in fl from the 
reference to the transformed lattice using an affine con- 
nection. This is done by minimizing the (positive) scalar 
quantity, 

n^fl i j 

x{rUO)-rimr (4) 

with respect to choices of affine e,y. Here the indices 
iandj = 1,2 (or x,y) and rn{t) and r^(0) are the i*'^ 
component of the position vector of the n*'* particle in 
the reference and transformed lattice, respectively. Any 
residual value of -DQ(r, t) is a measure of non-affineness. 

The rhombic lattice is a special case of the general 
oblique lattice - one of the five possible two dimensional 
Bravais lattices. In general, we need two order parame- 
ters (OP) to describe the transition between square and 
oblique lattices, i.e., between the space groups p4mm 
p2. These are the affine shear strain 63 — e.xy — ^yx and 
the deviatoric strain 62 = {exx~£yy)- Thus from symme- 
try considerations alone, we would expect to obtain four 
symmetry related product phases (25|. However, the mi- 
croscopic model used by us obtains a rhombic lattice for 
which 62 identically vanishes and the four equivalent vari- 
ants merge in pairs to give two symmetry related prod- 
ucts. To show this we have plotted in Figl3] the T — 
energy obtained for our model solid with the parameters 
a — 1, 113 — 0.2 and at p = 1.05 for various values of 62 
and 63 taking the square lattice as the reference. Apart 
from the minimum corresponding to the square lattice 
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FIG. 3: (color-online) Contour plot of the zero temperature 
energy per particle for p = N/V = 1.05, a = 1 and W3 = .204 
as a function of the OP strains (e2, 63) showing a metastable 
square minimum at (0, 0) and two degenerate, stable rhombic 
[oblique) minima at (0,±.18). 

ima is zero. It is therefore sufficient to use 63 as the sole 
order parameter (OP) for this transition [sj. 

Figure Hl^a) shows the nucleation and growth of 
the twinned martensite nucleus, following the lower 
temperature quench - we have plotted the best-fit 63 for 
snapshot configurations of = 110 x 110 particles at 
time steps of 2000, 3000, 4000 and 5000 Ai. The twinned 
structure of the nucleus composed of the two degenerate 
rhombi (characterized by positive and negative values 
of 63, separated by a sharp boundary) is evident even 
at the earliest time, and becomes more pronounced as 
time progresses. The constraint of fixed density forces 
a dynamical coupling between the affine OP strain and 
the affine non-order parameter (NOP) volume strain 
ei — + ej/j/j so that the transformation is also 
accompanied by a volume change, Fig|4l^b). As a result, 
as the transformation proceeds, more and more particles 
are pushed up against the surrounding untransformed 
square lattice which creates a jammed region at one 
end and an unjammed region at the other end of the 
anisotropic martensitic nucleus, Figdl^b). 

2. The dynamics of transformation is described by 
an affine OP strain (here, shear strain) characterizing 
the microstructure of the growing nucleus, and an affine 
NOP strain (here, volumetric strain), which is slaved to 
the former. 

We can now use the residual Df^, ^ to extract the 
spatio-temporal variation of any non-affine deformation 
that is produced during the transformation. To be able 
to distinguish between non-affineness arising from differ- 
ent components of the strain (shear or volumetric) distor- 
tion, we need to incorporate the notion of jamming in the 
definition of non-affineness In the context of gran- 
ular compaction[3^ and glassy materials[3^, jamming 
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FIG. 4: (color-online) Best-fit afline strains and residual non- 
affine deformations, obtained from MD simulations of parti- 
cles undergoing a square to rhombic transition at T = 0.1, 
for time slices (i) 2000, (ii) 3000, (iii) 4000 and (iv) 5000 At. 
Plots obtained by coarse-graining the A'' = 110 x 110 lattice 
to a 64 X 64 lattice, (a) Order parameter (shear) strain es, 
colors show 63 from —0.3 (black) to (brown) to 0.3 (yellow). 
The twinned microstructure is clearly visible, even at earlier 
times, (b) Non-order parameter (volumetric) strain ei, colors 
show ei from —0.5 (black) to 0.3 (yellow). The equilibrium 
value of ei is nonzero within the rhombic phase, in addition, 
ei appears at the two ends of the twinned microstructure due 
to elastic coupling to the order parameter es. (c) Non-affine 
deformation x- Colors show x ranging from —1 to 1. Note 
that X ^ at the centre of the growing nucleus at large times. 
The NAZs surround the growing nucleus and are created at 
and advected by the front. Jammed and unjammed NAZs 
occur at the 'top' and 'bottom' of the nucleus, respectively, 
sharing the same spatial symmetries as ei. 



this translates into computing the relative change in the 
distance between particles within U, in the direction of 
motion of the particles in the nucleus, denoted by A^; we 
may thus define a quantity x(r, t) = —D^ sign(A/), which 
takes both positive (jammed) and negative (unjammed) 
values. For the martensite nucleus, the jammed and un- 



vected by the transformation front. Note that the spatial 
symmetries of x are the same as that of ei at all times 
(compare FigslH^b) and (c)), and so we associate the 
non-affineness predominantly with the NOP or volumet- 
ric strain. This restriction of non-afhneness to the NOP 
strain alone, could be specific to the square-to-rhombus 
transition; in transitions between other structures, there 
could be a fair degree of plasticity associated with the 
OP strain too. We will return to this point in section IV. 
Consistent with geometrical theories, the NAZs are ab- 
sent at the twin interface; this interface is coherent and 
the twins are simply related to each other by an affine 
transformation. 

In contrast, Fig[5lja) shows the nucleation and growth 
of the ferrite nucleus, following the higher temperature 
quench - as above, we have plotted the best-fit 63 for 
snapshot configurations of = 110 x 110 particles at 
time steps of 8000, 10000, 13000 and 15000 M. The nu- 
cleus is composed of polycrystalline grains of the rhombic 
phase separated by large angle grain boundaries. As time 
progresses, the grains rotate with respect to each other, 
giving rise to large non-afhnc distortions even in the hulk 
of the nucleus. This is reflected in the large values of 
X in the bulk of the growing nucleus, Fig[5]^b). How- 
ever, a spatial average of the instantaneous values x and 
63 over a scale larger than the grain size, gives zero for 
both. Similarly, a time average of the local x S'lid 63 
over a window corresponding to typical grain reorganiza- 
tion times, gives zero for both. The plastic zone spreads 
throughout the product region causing extensive atomic 
rearrangements. 

We now take a closeup look at the NAZs - Fig[n)Ja),(b) 
shows snapshots of the atomic positions in the NAZs 
of the martensite and ferrite nucleus, respectively. It is 
clear that the atomic configurations in the NAZs are 
highly amorphous, with no clear relation to the reference 
parent lattice. It seems meaningless to describe NAZs 
in terms of a density of dislocations, since the reference 
state has no unique physical significance for character- 
izing the current state in the NAZs. Even if we were to 
describe the state of NAZs in terms of dislocations, the 
density of dislocations would be so high as to have over- 
lapping cores, thus rendering this language inadequate. 
It is more reasonable to describe the NAZs in terms of 
fluctuations in the local density 0(r, t) = {p{r,t)~po)/pQ, 
where p{r, t) — J2nen ^("^ ~ ^n{tY), and po is the average 
uniform density. Indeed in [6[, we had studied the 
dynamics of 0(r, t) in great detail and demonstrated 
its involvement with dynamics of transformation and 
microstructure selection. Here we find by explicit 
computation that (j) and x are related - localized regions 
with large </> correspond to large x and so on. 

3. Right from its initiation, the transformation 
is accompanied by non-affine deformations primarily 
associated with NOP (here volumetric) strain. The 
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FIG. 5: (color-online) Best-fit affine strains and residual 
non-affine deformations, obtained from MD simulations for 
a quench at T = 0.8, for time slices (i) 8000, (ii) 10000, (iii) 
13000 and (iv) 15000 /Si. Same coarse- graining as in Fig|4] 
(a) Order parameter (shear) strain 63, colors show 63 from 
—0.3 (black) to 0.3 (yellow). Note that unlike the marten- 
site, no clear spatial pattern in 63 can be discerned. The 
local structure of the product nucleus is polycrystalline, with 
individual grains which coarsen with time, (b) Non-affine 
deformation Xy colors show x ranging from —1 to 1. Note 
that non-affine regions are present throughout the interior of 
the nucleus, signifying extensive plastic deformation during 
growth. 



(b) 




FIG. 6: (color-online) Close up of a region from the nucleus 
(a) in martensite : corresponding to Fig|3][iii) and (b) ferrite : 
corresponding to FiglSjiv). The color code is as follows: black 
dots - untransformed regions; yellow and black circles - affine 
regions with -|-ve and — ve 63; red and blue circles - jammed 
and un-jammed non-affine regions. In (a) a similar non-affine 
region arises at the other (bottom) end of the twinned region 
(not shown). 



We will now show that NAZs are produced when the 
local volumetric stress exceeds a threshold value. We 
compute the instantaneous local stress from our MD sim- 
ulations by spatially averaging the generalized virial, 



F,;rr 



considerations of proper averaging and obtaining infor- 
mation over a fine enough length scale. We have chosen 
M = 100 as a compromise between these considerations. 
Further, in order to obtain good statistics for any time 
t, we average over many independent quench runs. Thus 
the spatio-temporal resolution of the computed a is not 
as high as the one for the coarse-grained strain e. Details 
pertaining to the evaluation of this quantity for our po- 
tential including three-body terms is given in Appendix 
A. 

We can now compute the local volumetric stress cti, 
affine volumetric strain ei and non-affine x, averaged 
over the coarse-grained cell r^M, at different times 
following the quench. This is plotted in Figl71[a),(b), 
where we have expressed the local stress as a fractional 
difference about the value of cti for ei = 0, viz., the 
undistorted region. The <Ji-ei plot shows a linear elastic 
regime for those coarse-grained cells where the strain 
ei is small; concomitantly the non-affine x is zero 
(Figl7Jc),(d)). Coarse-grained cells where ei is larger 
than a threshold, show yielding (nonlinear and erratic 
(Ti-ei) and appreciable plastic flow, x 7^ 0- We have 
verified that these coarse-grained cells showing plastic 
deformation are indeed the NAZs reported above. We 
now focus on one coarse-grained cell, and study the 
time development of ci, ei and x as the transformation 
proceeds (FiglS]). We find that at earlier times, the 
strains are small and the stress-strain response is elastic. 
Beyond a yield stress ctic, the stress-strain relation is 
nonlinear, giving rise to non-affine deformations X 7^ 0- 
Following yielding, the local stress eventually decreases, 
often exhibiting oscillatory behavior. We find that the 
threshold stresses ctic, when expressed as a fraction 
of the ambient stress is only weakly dependent on 
temperature. 

4. Non-affine deformations are produced when the 
local stress crosses a threshold. The threshold stress is 
only weakly dependent on temperature. 

In the next section, we will use the four principles high- 
lighted above to construct an elastoplastic theory for the 
dynamics of solid state transformations. These four prin- 
ciples are generic, and should not depend on the choice 
of potential, or the nature of the transformation. We 
will show that the development of the microstructure is 
crucially influenced by the dynamics of the NOP plas- 
tic strain associated with NAZs. In what follows we will 
develop this theory first in a general setting and then 
specialize to the particular case of the square to rhombic 
transformation studied in this section. 



III. ELASTOPLASTIC THEORY OF 
NUCLEATION DYNAMICS OF SOLIDS 
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FIG. 7: Local stress expressed as a fractional difference from 
the volumetric stress at ei = (Aai/(Ti(0)), plotted against 
the local strain at different times (symbols) obtained for 
the (a) quench at T = 0.1, averaged over 40 independent 
quenches, and (b) quench at T — 0.8, averaged over 35 in- 
dependent quenches, (c) and (d) corresponding plots of x 
vs ei. The linear Hooke-an regime represents the local elas- 
tic response at small local stress. Beyond a threshold, the 
system yields locally, giving rise to a nonlinear stress-strain 
behavior and simultaneously non-afRne deformations x > 
(red symbols). The regions in real space associated with this 
local plastic regime are identical to the NAZs. 
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FIG. 8: Time dependence of local stress A(Ji/(7i(0), strain ei 
and X near the growing nucleus after a quench at T = 0.1. 
Note that initially x = and the stress (apart from large 
statistical fluctuations) is proportional to strain. When the 
local stress exceeds a threshold, x begins to increase, and the 
stress versus strain is highly nonlinear. 



plastic strain to the free energy. 

The driving force for is the chemical potential gradi- 
ent 6J- /Se^ to form the rhombic phase and the dynamical 
equations take the general form[14]. 



*(eT,eT,eT;ey) = 



(6) 



The dynamics of the afline NOP strain is slaved to 
the OP strain. This takes the form of a local force 
balanceWOl, 



tensor as e = {ei^,ev}, where the affine OP or transfor- 
mation strain connects the parent and product lat- 
tices, and the NOP strain ey is split into an affine part 
and a non-afRne or plastic strain, ey — ey ~ ty (ey can 
have many components, we have dropped the tensor in- 
dices for clarity). Note that our association of plasticity 
with the NOP strain alone, follows from the previous 
section, however as we remarked, this could be special to 
the square-to-rhombus transformation. In general, there 
could be plastic deformations associated with the OP 
strain too; we will comment on this in section IV. 

The transformation is described by a free-energy func- 
tional 



V • cry = 



T = j FT{{e^})+a{v4f+p{Ve^f + ^\ke^~. 

(5) 

with FT{{e^}) having three minima corresponding to the 
parent phase = and (symmetry related) variants of 
the product phase, ^ 0. This choice of T is dictated 
by simplicity, more complicated forms having symmetry 
allowed nonlinear cross-couplings can be envisaged, but 
these will not change our framework. In the absence of 
plastic deformations, a variation of (O with respect to 
Cy, gives the desired relation connecting the OP and the 



(7) 



where the local NOP stress cry is related to the instan- 
taneous equilibrium value of the NOP strain. 



cry = 



d4 



(8) 



The dynamical equations for the plastic NOP strain ey 
are constructed phenomenologically following the princi- 
ples listed above. We include the physics of threshold 
stress and yield flow, by a constitutive relation between 
stress and strain rate. 



■ p 



1 , 
h 




cTo)" if g{av;<7Vc) > 
otherwise (9) 



where h, 9, ao and avc are material parameters and 
g{(7v',crvc) is the appropriate (material dependent) 
threshold or yield criterion [2ll [22 |. which in principle 
can incorporate history dependence. In the next sec- 
tion, we will make the simplified 'Newtonian' ansatz, 
9 = l,aQ — 0, for this case the parameter h is related 
to the relaxation time of plastic flow. The ratio of the 
stress to strain rate, a/iy = rj, is a bulk viscosity. At 
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locally decrease and increase once the solid yields, giving 
rise to oscillatory behavior, as seen in FigO The form 
of the constitutive relation [9] precludes the possibility of 
"creep" since ey vanishes at zero stress. This motivated 
by our MD results as shown in Fig[Sl 

In addition, owing to plasticity, the local afFine strains 
do not satisfy the usual St. Venant's compatibility[3| 
— instead, the amount of incompatibility is exactly ac- 
counted for by the amount of plasticity generated. This 
implies that the local St. Venant's condition should be 
rewritten as [4311 . 



V X (V X e)^ 







(10) 



where e is the total strain, which includes ey. In regions 
where the local plastic deformation is zero, this reduces 
to the usual St. Venant's compatibility condition. 

The equations of constraint, ([7]) and (fTU|. are used to 
express ey in terms of e^; equations, ([6]) and ([9]), are 
then used with appropriate initial conditions to describe 
the elastoplastic theory for the dynamics of solid state 
transformations. 

As we highlighted in the previous section, the nucle- 
ation process is heterogeneous. This implies that since 
equations ^ and ^ are deterministic, one needs to in- 
troduce 'seeds' in order to initiate nucleation events. In 
our elastoplastic theory, this can be introduced as initial 
conditions either in the local strain or the local stress 
field cry. A random (initial) distribution of the stress 
field cry, can produce local plastic strain Cy, if the lo- 
cal cry seed is larger than the yield stress. This will in 
turn nucleate the transformed solid, via e^. The results 
obtained from solving the dynamical equations must be 
then averaged over realizations of the quenched random 
stress field. A simpler strategy is to directly introduce a 
seed in the transformation strain, , by creating a small 
twinned region with a single twin boundary (since the 
dynamics of is conserved). This allows us to follow 
the subsequent dynamics in precise detail and does not 
require any averaging over noise realizations. 

We use the elastoplastic model to address two separate 
but related issues. We will first determine the late time 
morphology and micro structure of the growing nucleus 
following a quench across the structural phase boundary 
and construct a dynamical phase diagram akin to Fig[2] 
of our MD simulation, or Fig. 4 of Q . We compute the 
shape of the growing nucleus from the shape asphericity 
A = {Xi — A2)/(Ai -I- A2), where A; are the eigenvalues 
of the moment of inertia tensor of the nucleus [27|. By 
quenching into different regions of the dynamical phase 
diagram, we will study the dynamics by which specific 
microstructure gets selected within our elastoplastic the- 
ory. 

A brief comment, before we discuss our explicit com- 
putation : strain-only theories [14l| set Cy = 0' ^'^d use ^ 
and (fTU|) to eliminate ty in terms of e;^, leading to long- 
ranged interactions and dissipation in e^, both of which 



in the strain-only description; the size of the twins is set 
by elastic parameters alone^l^]. These theories cannot 
describe the occurrence of the ferrite. In our elastoplas- 
tic description, inclusion of local plastic deformation in 
the form of ey has two effects — it screens and isotropises 
the non-local interaction and dissipation kernel. This is 
ultimately responsible for the destruction of the twin pat- 
tern, resulting in a ferrite. Further, the size of the twins 
depends on elastic, as well as plastic parameters. We pro- 
vide a detailed analysis of these effects in a forthcoming 
publication. 



A. The square to rhombic transition 

We now present an explicit calculation for the particu- 
lar case of the square to rhombus transformation in two 
dimensions, studied in section II. As mentioned there, 
the square to rhombus transformation is a special case 
of the square (p4mm) to oblique (p2), and is thus de- 
scribed by two order parameter strains characterising the 
4-degenerate product phases. In general, one can con- 
struct a Landau theory for this transition(25l. [2^ using 
terms upto sixth order in the OP strains = {62,63} 
and quadratic in the NOP strain ey = ei. We then 
decompose the NOP strain into a slaved, affine NOP 



strain ey 

p — p 
e'y = ef 



A — „A 



Ci, and a dynamical non- affine NOP strain 
e{ , enabling the total NOP strain to be written 
as, ei = — ef . With this decomposition, we can then 
proceed with the general treatment outlined in the last 
subsection. 

However to compare with the results and phenomenol- 
ogy of the MD simulation of the model solid described by 
the microscopic potential ([2]), ([3]), it is more convenient to 
look at a restriction of this problem. Recall (see, Figl3]) 
that our microscopic potential supports two rather than 
the possible four product minima. This implies that for 
this choice of potential, there is therefore only one min- 
imum in the 62 direction; it thus suffices to retain upto 
quadratic terms in 62 in the strain free-energy functional. 
The minimal free-energy functional, sufficient to describe 
this square to rhombus transition is given by, 

^ = \j dxdy [ai(6i -fef)2 + aae^ + 036^ (11) 

+ Ci(V(6i + 6f ))2 + C2(V62)2 + (V63)2 - 6^ + 6^] 

in terms of the OP strains 63 and 62 and the NOP strain 
61. The only other term to quadratic order in 62, viz., 
63 e| has also been dropped since it does not influence 
the phase transition. Note that this form of the free- 
energy functional can be recast just as in ([5]). The three 
elastic constants 01,02,03 define the linear elasticity of 
the square phase. The coefficients of the quartic and 
sixth order terms as well as that of Ves can be scaled 
to unity by rescaling 61, J- and the spatial coordinates 




FIG. 9; Pictorial representation of the dynamics of the plastic 
strain ef as a function of the local volumetric stress ai given 
by @ for a typical flowing solid (solid line) and the simplified 
form (fTG]) used in our computations (dashed line) . 



crystal is metastable and the rhombic crystal is stable at 
equilibrium. 

The affine NOP strain is slaved to the OP strains; we 
make use of the two conditions ([7]) and (|10p . to express 
in terms of 62 and 63. Mechanical equilibrium implies 
V • CT = 0. The modified St. Venant condition (fTOt reads, 



4V,Vye3 = 0, (12) 



which includes both affine and plastic deformations. The 
relation between the total NOP strain ei and the OP 
strains is most conveniently expressed in k— space, 



,(k) = Q.3(k)e~3(k) 



(13) 



with the kernels, 



Qi3(k) 



4a2 



2^3 hxky 



ai + 02 k? 



= qi3 



(14a) 



and 



Qi2(k) 



as 



2a2 k'l 



kl 



2ai + a3 k"^ 



= -qi2 



k^ ky 



(14b) 



In effect, the above equations of constraint, connect the 
instantaneous to the dynamical ef , 62 and 63. 

We now have to specify the dynamics for the OP 
strains and the plastic NOP strain. At this stage we make 
the approximation of replacing the value of the OP 62 by 
its value at equilibrium, i.e., 62 = for all times. This 
greatly simplifies the calculation without changing the 
physics. The dynamical equation for the affine OP strain 
Q 63 may be derived from Newton's laws[l4| incorporat- 
ing dissipation via a Rayleigh dissipation functional [40|. 



FIG. 10: Results from the numerical solution of the dynamical 
equations in a 128 x 128 grid with a time step St — .002. The 
parameters for this calculation are ai — 100, a2 — 1,03 — 
.01, 7 = 5. (a) Shape asphericity A of the growing nucleus 
as a function of time t from the elastoplastic model - bold 
line martensite and thin line ferrite. (b) Dynamical phase 
diagram in the aic — hi plane (expressed in units of aiAe^ 
and 7, respectively), starting from the same (elliptical) initial 
seed. The shape asphericity A at later times and the order 
parameter strain ea have been used to determine the phases. 



where 7 is a solid shear viscosity. 

The dynamics of the plastic NOP strain are determined 

by m, 



1 

hi 

if \ai \ > aic 
= otherwise 



(16) 



where we have, for simplicity, chosen 9 = 1 and a simple 
threshold criterion, with yield stress aic- 

We have solved the dynamical equations (|15p and (fTB|) 
by discretizing in space and time. The initial conditions 
are chosen from an appropriately seeded (metastable) 
square phase (see discussion at the end of previous sub- 
section). Details of the numerical computation are given 
in Appendix B. Note all quantities have been made di- 
mensionless. 

As in the MD simulation, we first obtain a dynamical 
phase diagram demarcating the regions where a marten- 
site or ferrite is obtained upon quenching. This is done 
by studying the shape of the nucleus, in terms of the 
shape asphericity A, and twinning of the microstructure, 
in terms of the affine OP strain 63. We focus on a special 
cut in parameter space; we fix the coefficients appear- 
ing in (jl2p and 7, and explore the dynamical phase dia- 
gram in the plasticity variables, aic-hi plane (Fig llOr b)). 
The threshold stress uic is expressed in units of oiAe^ 
(the affine stress at the structural transition), and hi in 
units of P (where ^, thickness of the twin interface, 
has been taken to be 1). In FigfTUTa). we plot A as a 
function of time for an anisotropic (twinned) nucleus and 
an isotropic (untwinned) nucleus; the value of A at late 
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containing the martensite and ferrite, FigfTUrbV 

We discuss several novel features of this dynamical 
phase diagram. For instance, even when the threshold 
stress (Tic is zero, a martensite can form if the plastic- 
ity relaxation rate is small (large hi) compared to the 
rate of growth of the nucleus. This feature was already 
present in our earlier calculation 6] , where the dynamics 
of the local density fluctuations determined the selection 
of microstructure, and is an inescapable feature of real 
martensites[4^. The phase diagram, Fig fTUf b). is con- 



structed for a fixed value of under-cooling. As the degree 
of under-cooling changes, the phase boundary changes 
slightly, but not a whole lot. More significantly, the plas- 
ticity relaxation time hi increases with the lowering of 
temperature. Thus by starting out in the ferrite phase, 
one can cross the phase boundary into the martensite 
by simply lowering the temperature, identified as the 
martensitc-start or Ms temperature. In addition, there 
is a well defined plateau yield stress over three decades 
in hi, suggesting that the yield stress is independent of 
temperature over this range. This is consistent with our 
MD simulations. Finally, it must be noted that the dy- 
namical phase diagram FigfTUTb) is constructed from the 
nature of the first critical nucleus that forms. In a macro- 
scopic sample, a ferrite nucleus may eventually nucleate 
and grow even in the martensite phase, once the plas- 
tic strain ef gets enough time to relax. This is consis- 
tent with our MD simulations and agrees with results of 
isothermal quenching experiments in real materials @. 

Having displayed the dynamical phase diagram we can 
perform quenches to the martensite and ferrite phase, 
and study the time development of the profiles of 63, cti 
and ef (Figs. [Til [13)) . Using plasticity parameters cor- 
responding to the martensite phase, (Fig. [TT|) shows the 
temporal evolution of a twinned nucleus, in perfect anal- 
ogy with our MD simulations. The nucleus initially grows 
parallel to the twin boundary fFig. [TlT i)). while the stress 
CTi approaches the threshold at the growing tips. As a re- 
sult the plastic strain ef gets to be large at these tips. 
As the nucleus grows, these highly stressed and plastic 
regions are advected by the growing tips. In the interior 
transformed region, the stress relaxes to a sub-threshold 
value, and the plastic deformation goes to zero. The se- 
quence of events exactly mimic the dynamics of the ey, 
X and NAZs of Fig. [H To study the time evolution of 
the NAZs in more detail, we focus on a single cell, fl, 
midway within our computation box, along (but not on) 
the twin boundary of the initial seed within the untrans- 
formed square lattice. With time, the growing tip of the 
nucleus approaches, and then sweeps by f2, in the process 
transforming it into the triangular phase. This situation 
is analogous to that shown in Fig[8lfrom our MD simula- 
tions. We plot the local ei,ai and ef at 57 as a function 
of time t in Fig[T21 As in Fig[51 initially fl, which lies 
ahead of the approaching transformation front, begins to 
deform elastically due to stress generated at the growing 




FIG. 11: (color-online) Time development of the afiine strain, 
stress and non-afRne strain following a quench into the 
martensite phase for 256 x 256 cells, at (i) t — 40, (ii) 
t = 350 and (iii) t = 800, starting from an initial elliptical 
nucleus with a single twin boundary. The plasticity param- 
eters |(Tic| — 1 and hi = 1 while the rest of the parameters 
are as before, (a) Profile of afBne OP strain 63, showing the 
initial growth parallel to the twin interface, followed by the 
dynamical addition of twins. Colors: yellow to black maps the 
range — 1. < 63 < 1. Brown region denotes retained austenite, 
63 = 0. (b) Corresponding profile of the local stress ai. The 
local stress is concentrated at the tips of the growing front 
where it approaches the threshold value aic (i) and (ii). In 
the interior of the growing nucleus cri relaxes to zero. Sub- 
sequently in (iii), the local stress gets large in regions where 
the new twins are being accommodated. Note the variation 
in the signs of the stress in the direction along which new 
twins are added. Colors: yellow to black maps the range 
— 1. < (Ji < 1. (c) Corresponding non-afhne strain ef , show- 
ing the initial advection by the transformation front, and its 
dynamical emergence as subsequent twins are added, ef ap- 
pears in regions where ai ~ uic and cri ^ — o"ic. Colors: 
yellow to black maps the range —.01 < ef < .01. 



approaches H, ai rises and tends to cross the threshold, 
(Tic. At this instant, ef begins to form reducing ai to a 
value below aic- As the nucleus grows further, cri within 
increases again - and the process repeats producing a 
local stress which oscillates rapidly in time. These oscil- 
lations result from cooperative jamming and unjamming 
events caused by alternating build-up of ai due to inter- 
face motion and its relaxation by creation of NAZs 3^. 
The region of high (and oscillating) local stress and the 
NAZ travels with the growing tip, being advected by the 
moving transformation front. The strong resemblance 
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FIG. 12: Evolution of the local total strain ei, the local 
stress (Ji and the non-affine part of the strain ef for a cell 
n — (128, 100) on the twin axis following the quench into 
the martensite phase shown in Fig llll Note the initial linear 
regime, when ai cx ei and ef = 0, followed by oscillations in 
CTl and the creation of ef as the stress cti rises to the thresh- 
old value CTic = 1. The resemblance with Fig|8]from our MD 
simulations is quite apparent. We have multiplied the strains 
by 100 in order to plot them to the same scale. 



the bulk of the product phase. 

The subsequent dynamics, Fig lllf iii). goes beyond the 
time scales accessed in the MD simulation. The dy- 
namics now proceeds perpendicular to the twin interface, 
adding new twins (symmetrically disposed) as time pro- 
ceeds. The production of new twins with a fixed width, is 
a consequence of the anisotropic, non-local interactions 
connecting spatially separated regions with nonzero 63. 

Note that the affine NOP strain mediates the 
anisotropic, long-ranged interactions[14]; the presence of 
the plastic NOP strain screens this interaction, making 
it short-range (Appendix B). The emergence of the OP 
strain in the form of twins, leads to an increase in the 
local stress ai, which in turn generates plastic flow on 
crossing the threshold (Fig fTTT c)). The plastic strain, 
once produced, reduces the value of the total NOP strain 
and, therefore, that of the non-local interaction. In the 
case of the martensite, this reduction is not complete. 
The stress cti decreases to zero in the interior of the 
martensitic nucleus, so that any given region undergoes 
the same sequence of transformations : untransformed 
— > elastic distortion non-affine —^ transformed. This 
temporal sequence is also seen in our MD simulations. 

We now use plasticity parameters corresponding to the 
ferrite phase; FiglT2] depicts the time evolution of an 
isotropic, polycrystalline ferrite nucleus, starting from 
the same initial conditions as above. For small values 
of (Tic, the plastic strain ef is produced readily, and on 
an average largely cancels out the effect of the affine 
NOP strain. This significantly reduces the magnitude 
and range of the non-local interactions, which were re- 
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FIG. 13: (color-online) Time development of the affine strain, 
stress and non-afSne strain following a quench into the ferrite 
phase, at (i) t = 40, (ii) t — 300 and (iii) t — 1600, starting 
from an initial elliptical nucleus with a single twin boundary. 
The plasticity parameters ctic = and hi = .1 while the 
rest of the parameters are as before, (a) Profile of affine OP 
strain 63, which shows the initial elliptical nucleus growing 
approximately isotropically. Colors: yellow to black maps the 
range — 1. < 63 < 1. (b) The local stress cti rapidly relaxes to 
zero, in this case, with only a small residual value remaining 
in the vicinity of the original seed. Colors:yellow to black 
maps the range —1. < cti < 1. (c) Corresponding non-affine 
strain ef , showing its invasion into the 'bulk' of the growing 
polycrystalline nucleus. Colors: yellow to black maps the 
range —.1 < ef < .1. 



the interior of the growing nucleus, which leads to an 
invasion of the non-affine strain ef into the 'bulk' of the 
isotropically growing, polycrystalline nucleus. This re- 
sults in incoherent grain boundaries in the interior of the 
growing nucleus. Note that the symmetry of the pattern, 
FigHni arises because the evolution equations are de- 
terministic; any noise would destroy this symmetry and 
make the grain boundaries rough and orient randomly. 
The sequence of events then exactly mimic the dynamics 
of the e^, X and NAZs of Fig. [3 

We believe we have successfully constructed a general 
elastoplastic description for the dynamics of solid state 
transformations, which is capable of describing different 
microstructures and addresses the issues of microstruc- 
ture selection. The qualitative picture that emerges from 
the elastoplastic model closely resembles our MD simu- 
lation results. In the next section, we will present some 
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IV. DISCUSSION 

By constructing a theory for microstructure selection, 
we successfully bridge two apparently disparate descrip- 
tions of the dynamics of ferrites and martensites, within 
a unified framework. This has been achieved at the cost 
of enlarging the space of dynamical variables to include 
non-affine deformations. Have we lost some of the special 
features of martensites in the process ? What implica- 
tions does our description of solid state nucleation and 
microstructure selection have for conventional nucleation 
theory? 

1. Emergence of average compatibility from the elasto- 
plastic dynamics : Our elastoplastic theory pro- 
vides an understanding of how the strict local com- 
patibility of strain-only theories [l4|, [l6| can be rec- 
onciled with the average compatibility of geomet- 
rical theories [l3|, when describing the dynamics of 
the martensite. We find that in the regime where 
the martensite obtains, our dynamical equations 
reduce to the equations of the strain-only theo- 
ries, when we coarse-grain over a scale A corre- 
sponding to the size of the NAZs. Simultaneously, 
the equations of constraint, viz., the modified St. 
Venant's condition (fTU|) . reduces to the usual 
St. Venant's elastic compatibility, provided we 
coarse- grain over the same scale A. This can be 
seen by explicitly writing out p^ . 

V'ef - (V^ - V2)e2 - 4V,Vye3 = V^ef . 

In Figim we have re-plotted ef for time t — 800 
(FiglTIJc)(iii)) in order to show the NAZs in detail. 
Remarkably, the martensite plates are accompanied 
by patches where ef and alternate in sign. If 
the system is coarse-grained over a distance A the 
effect of these patches cancel and we recover the 
usual elastic compatibility condition. This emer- 
gence of compatibility upon coarse-graining over 
the scale of the plastic zone was noted in [H . As the 
critical stress ctic increases, the size of the NAZs 
decreases and so does the coarse-graining length 
scale A. Note that the coarse-graining appropriate 
for the emergence of the usual elastic compatibil- 
ity, does not wash out the twinned microstructure 
of the martensite; the 'phase' of the averaging is so 
as to produce a planar martensite-austenite inter- 
face as shown by the dashed line in Fig[TJ 

2. Reversibility of martensitic transformations : How 
do we reconcile the inevitable creation and evolu- 
tion of plasticity in the form of NAZs with the ap- 
parent reversibility of martensitic transformations, 
as observed in shape-memory alloys ? Here we will 
provide some preliminary comments, which will be 
taken up in greater detail later[4i|. In essence, mi- 
crostructural reversibility in martensites, is related 




FIG. 14: (color-online) The plot of the plastic strain ef for 
a growing martensite at t = 800 (same as in Fig lllf cKiii')'). 
Note the presence of NAZs with alternating signs of ef accom- 
panying the growing martensite phase - see region on top the 
within the dashed lines, a similar region (not marked) exists 
near the bottom of the figure. This figure needs to be com- 
pared with Fig[T] When the strain fields are coarse-grained 
over the length scale A the contribution from these alternating 
patches cancel, and full elastic compatibility is restored. 

NAZs, is that it is largely associated with the NOP 
sector, which in turn is slaved to the transformation 
strain. Indeed even the slightest amount of plastic- 
ity in the OP sector, would make the transforma- 
tion irreversible. This is apparent in martensites 
involvin g F e alloys, which do not exhibit shape- 
memory [42|. Within our own model system, a deep 
quench to the a — 0, V3 = region produces a 
triangular solid which is not related to the par- 
ent square lattice by a group-subgroup relation[25|. 
During the reverse transformation, therefore, there 
is no unique parent lattice that the system can re- 
vert to. This produces non-afHneness in the OP 
sector due to a multiplicity of affine paths and de- 
stroys reversibility. Our results related to such non- 
reversible transformations will be published else- 
where. 

The other relevant feature exhibited by the NAZs 
associated with martensites, is the special nature 
of the particle trajectories ja, [13 • Particles in the 
NAZs formed during martensitic growth, move bal- 
listically and in a coordinated manner. It is these 
two properties of the NAZs discussed here that ul- 
timately renders the square to rhombus marten- 
sitic transformation reversible, in spite of signifi- 
cant transient and localized plastic deformation. 

3. Inconsistency with Ostwald's step rule : Ostwald's 
step rule of 1897, states that "the phase that nu- 
cleates need not be the stable phase, but the one 
that is closest in free energy to the parent phase 
. . . ". This rule has been interpreted by Stran- 
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ated. While it is easy to appreciate the apphca- 
bility of this rule for phase transformations in sim- 
ple systems having a uniquely defined barrier cross- 
ing event along the path of the transformation, it 
is more difficult to apply such considerations to, 
say, atomic rearrangements and the generation of 
NAZs where many barriers with different attempt 
frequencies may be involved j46l|. As discussed in 
section III, the selection of microstructure depends 
both on parameters in the free-energy functional 

and dynamical parameters in (jTS]) and 
This is explicitly shown in Fig. 1101 where the mi- 
crostructure depends on the plasticity dynamics, in 
terms of the yield stress aye and viscosity h. 

Randomness and heterogeneous nucleation at defect 
sites : The results presented in section III A, were 
obtained with initial conditions corresponding to a 
small elliptical nucleus and the choice of dynam- 
ical parameters corresponding to martensitic and 
ferritic growth. An alternate initial condition for 
the nucleation dynamics is to prescribe a spatially 
random stress profile, e.g., a random a[. In real 
materials, this would correspond to frozen in de- 
fect structures. This quenched random stress would 
add to the internal stress so that the total stress 
CTi = aie° -t- a[. The intial stages of the dynamics 
of nucleation and growth is sensitive to the initial 
distribution, which we take from a Gaussian distri- 
bution with zero mean and width 5*1. 

For instance, if 5*1 <C CTic, the threshold, inhomo- 
geneous elastic strains in ei develop which initi- 
ate nucleation of the product, via its coupling to 
63. As the dynamics proceeds, local (total) cti 
gets enhanced at the transformation fronts, giving 
rise to plastic deformations. On the other hand, if 
Si > uic, the threshold, then local regions can de- 
velop appreciable plastic strains, Eqn ll6l The elas- 
tic coupling between ef and 63 through the func- 
tional derivative in Eqn [13 then causes the nucle- 
ation of the product phase (see Appendix B). The 
lag time for nucleation is appreciably faster com- 
pared to when Si <C ctic- The qualitative features 
of the phase diagram, FigUHl are unaltered by the 
presence of quenched random stress fields, though 
it 'enlarges' the regime over which ferrite phase ob- 
tains. The ferrite microstructure is altered; the 
grain sizes shrink and the grain boundaries thicken 
with increased randomness. The martensite mi- 
crostructure, at the scale of the twin pattern is how- 
ever unaltered, though again the the distribution of 
martensitic grains gets smaller with increased dis- 
order. This robustness of martensitic patterning 
over the scale of the twins is significant, and we 
wish to revisit this aspect in a detailed study. 

Dynamical phase diagram and TTT curves : 




FIG. 15: (color-online) Plot of the OP strain 63 for time t = 
1000 for a Martensite (a) with aic/aiAei = 2500 and for a 
ferrite (b) with aic — 0; /11/7 = .1 for (a) and (b). The 
color scheme is as in Fig lllf a'). In both these calculation, 
a quenched, Gaussian, random stress ai(r,t) with zero mean 
and variance Si = 10 was used to nucleate the product phase. 
Note that quenched randomness does not affect the overall 
characteristics of the Martensite. The ferrite grains however 
become much smaller with grain boundaries which are broad 
compared to the size of the grains. 



Transformation (TTT) diagram, constructed in the 
form of contour plots of the proportion of each con- 
stituent (martensite or ferrite) as a function of time 
during an isothermal transformation at different 
quench temperatures [l|, Q. Figure [TUl is the dy- 
namical phase diagram computed within our elasto- 
plastic theory; this can be converted to the typical 
TTT curve, provided we know the temperature de- 
pendence of aic and hi. Such a temperature de- 
pendence may be put in phenomenologically, as in 
Fig. 4 of 0, or obtained from a first principles 
non-affine field theory. 

6. Future work : It should be possible to extend our 
elastoplastic theory to include the effects of im- 
purities such as interstitial carbon in Fe (as in 
steel) [13], which undergoes significant non-affine 
deformation Q. Interstitial carbon represented by 
a diffusive concentration field ^, would enter into 
both the dynamical equations for the affine and 
non-affine strains. Thus any attempt to understand 
microstructure selection in systems such as steel, 
would involve the study of the coupled dynamics of 
the affine strain, non-affine strain and concentra- 
tion field ip. 

Finally, apart from restoring the full tensorial charac- 
ter of the elastoplastic description, we need to explore in 
greater detail, the consequences of general thresholding 
and yield flow in the plasticity dynamics accompanying 
solid state transformations. In a later paper (41i|. we will 
discuss the dynamical response of the transforming solid 
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APPENDIX A 

In this appendix, we reproduce the formulae used 
to calculate the energy, forces and stresses for the 
three-body potential used. The form of the three-body 
potential ensures that these quan tities are computable 
using pairwise functions alone [33j. This greatly speeds 
up the computation. 

Energy : The three-body part of the energy is given by 
^3 = ^ ^ i'3{rij,rjk,rkz) (Al) 



2^4 



^ 2 (sin^(e'jifc) cos^(0jife) 
4sin^(%fe) co^'^{djik)) fijfik- 



Now define 



Xij / 



and 



so that 



Using the above definitions and the quantities 

= x^jij^jfij, 
= x^jy^j{x^j — yj^j) fij, 
9ij{^) ~ x^jVij{xfj—ylj)fij, 

we get E-i = V3 Y,i Si with 

- 16 [G,(3)F, + 32G,(3)' + 2G,(4)2 
+ G,(l)2 - 16Gi(3)G,(l) 

- 4G,(5)G,(2) + 16G,(5)2] 

where Gi{n) = J^j^i 9v{n) and F^ = J^j^i /y - 

Force: The three-body forces can be found by taking 



rij and evaluating the quantities 



dX'i'i 



%j-(2) 



dxij 



%j(4) 

dXij 



= 2 [iyy2-(l -i^^ ) 

~ ^ijVij] fij/'''ij 

= [i^x^yij - y^j){l - ijj) 

— Xijyij{x^j — 3ij^yf^)] fij /rij 

+ 2i-jyjj(i-j - y1j){nj - ro). 



+ 24y^,(r„- 



[{Qx\jy% + 2i.,jy% 



- %~x%y%){X-~xl)-{^y%~xl 

+ 2yijXij — 8yijXij)xijyij~^ fij/^i, 

— {Sx^jijij — 5Xijyfj)\ fij /rij 
+ 2itjy|,-(i?j - yfj)inj - ro). 



the force acting on the particle 'i' in the x-direction owing 
to the three-body interaction can be written as 



%j(5) 

OXij 



(A2) 



where 



H^j = 4U3 



Ox in 



^ fik + 2xij {rij - ro) ^ gifc(l) 



dxi 



E^^-(l)-2^E5-(2) 



k^i 



k^i 



- 16U, 



%j(3) 
dx 



64 



dxi 



V y^/ ^ ^ 2ij^- (r,j - ro) ^ gik (3) 

k^i k^i 
Y.9^k{S)+A- 



k^i 



dx 



k^i 



+ 2«ii5:,..(i)-i6 



- 4 



dxi 

%j(5) 
dxi 



k^i 



dxi 



k^i 



k^i 



^ %.(2) 

3 Xij 



k^i 
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The j/-component of the force 



(A3) 



where Hfj is evaluated using expressions similar to that 
for H!^j given above. 

Stress : The contribution of the three-body interaction 
to the virial stress can now be calculated using the force 
components, 



(A4) 



where a, — x, y. 



APPENDIX B 



In this appendix we give details of the numerical solu- 
tion of Uni) and used to obtain the results of section 
III. We have used a simple real space scheme for the so- 
lution, discretizing the partial differential equations over 
a lattice of square cells of size 6x = 1. We have used 
128 X 128 cells to obtain the phase diagram fFig llOp and 
256 X 256 cells for FigslUllJ and [13l The initial value 
problem in time is solved using an Eulcr scheme with a 
time step of St = .002 which is sufficient to avoid nu- 
merical instabilities. Below we give the sequence of steps 
involved in the iteration of the discretized equations. 

Step 1. We start with initial values for 63, its time deriva- 
tive 63 and ef defined over all the cells in our lattice and 
for time t. First, we need to compute the (slaved) affine 
strain from the OP strain 63 by solving p^ ). The 
affine strain ef together with the known ef determines 
the total strain ei at t. In real space, [Hb becomes. 



^2 

ei{x, y) = 9i3^^e3(x, y) 



(Bl) 



Equation IBll is the Poisson equation for the charge den- 
sity, 

Pi{.x,y,{ez}) ^ qiz^-^e:i{x,y) (B2) 

and we need to solve it for the (Dirichlet) boundary con- 
dition ei — *■ for x,?/ — > 00. This is done by discretiza- 
tion in real space and by using an iterative scheme with 
a small over-relaxation and a convergence criterion of 1 
in 10^(4^. For convenience in what follows, we refer to 
this solution using the notation, ei = 'P{{pz}). 

Our numerics can be checked for accuracy by compar- 
ing the results with that of simple choices for 63 for which 
Ci (same as ei) may be obtained analytically. For exam- 
ple for 63 which is nonzero only within a square of size 
2a, viz., 



{Q{x) is the Heaviside step function) is the 

electrostatic potential for a set of four charges 
+eo,— eo,-feo and — eo at the vertices of a square 
(a, a), (a, —a), (—a, —a), (—a, a) in two dimensions. This 
is given by. 



ei{x,y) = 



£0 
2 



In 



[{x - af + {y~ aY 
[(x-ha)2 + {y- ay 

a^ + iy + aY 



[{x - a)2 + (y + a)2 



(B4) 



On the interfaces y — ±a and x — ±a it is easy to see 



that 



zbx (and ±y respectively) except near the cor- 



ners where there are weak logarithmic singularities. In- 
cidentally, this linear approximation for is the same 
as the initial value of the density fluctuation (j){r, 0) used 
in Refs.[3 and 0]. The presence of ef modifies the in- 
terfacial energy of a rectangular nucleus of martensite of 
length L and width W containing N twins leading to the 
experimentally observed scaling law L/N ^ W^^"^^. 

Our numerical result for for the choice of 63 given 
in (|B3p reproduces the analytic form (|B4[) to within a 
few percent. 

Step 2. Knowing the strain ei at time t, we next update 
the plastic strain to the next time step t + St hy iter- 
ating For this, the local stress ui = aici (plus any 
external stress if present) is obtained for all the cells and 
is then used as input to pB]) . 

Step 3. Lastly, we have to update 63 and 63 for which 
one needs to compute the functional derivative. 



ST _ ST3 ^ SJ^i 
Se^ Sez Sez 



(B5) 



where the two terms on the right hand side represent 
functional derivatives of the parts of the free energy ([T2|) 
involving only 63 and ei respectively. The first term is 
straight forward and is given by, 

^ = -V^es + e3 - 4(63)=^ + 6(63)^ (B6) 

and the second term, after some algebra can be shown to 
be. 



5_Ti 

Sez 



where, 
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(B7a) 



(B7b) 



Note that (|B7p involves the total ei which includes both 
affine as well as the non-affine strain. Even if 63 = to 
begin with, at subsequent times 63 may be created due 
to the presence of nonzero ef . One encounters such a 



Pi{x,y,{ez}) = 9i3^^ei(a;,y) 
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is similar to the seeding of the martensite nucleus with 
a single point-vacancy as in our MD simulations. Sec- 
ondly, the resulting form for the functional derivatives 
are highly non-local since they involve repeated solutions 



of the Poisson equation. However, if ef is large the to- 
tal NOP strain ei vanishes and the non-local coupling 
between spatially separated regions of the OP 63 disap- 
pears. 
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